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The implications of the most recent experimental results in the B and K meson systems on 
the CKM paradigm are investigated by means of a global fit to the theoretical predictions 
of the Standard Model. We advocate the frequentist approach 7?fit, which is implemented in 
^ I the CKMfitter package. Within this approach, constraints on the CKM parameters and other 

OO ■ quantities of interest are obtained, and extensions of the Standard Model are investigated. 

OO ■ LAL 02-79 

o . 

^ , 1 Introduction 

o 

c~| ' The electroweak sector of the Standard Model (SM) has been intensively studied during the last 

decades, in particular at LEP and SLD. CP violation in flavour changing processes is the next 
Q^l sector to be investigated. CP violation was first observed in the Kaon systemlil, and recently 

established in the i?-meson systemH, but within the SM it cannot produce a baryon asymmetry 
as large as what is observed in the universe a. A variety of models beyond the SM introduce 
new sources of CP violation: therefore, CP violation is an interesting place to look for new 
^ ' physics, and particularly in the i3-meson system where there are numerous large CP- violating 

observables. 



Within the Standard Model, CP violation is accounted for by one phase ^ in the CKM 
matrix&H which describes the charged current couplings among left-handed quark fields. With 
three generations, the CKM matrix is parameterized by four parameters (including the CP 
violating phase). We use the Wolfenstein parameterization 0, which, following the observed 
hierarchy between the matrix elements, expands the CKM matrix in terms of the four parameters 
p, r], A and A (A ~ 0.22 being the expansion parameter). Presently, A is known at the 1% level, A 
at the 5% level, whereas p and r] are the least known CKM parameters. Many relations between 
the matrix elements follow from the unitarity of the CKM matrix. In particular, the unitarity 
relation linking the first and third rows of the matrix is usually displayed as the Unitarity 
Triangle (UT) in the complex plane {p,f]) ^ 

Investigating the charged current couplings among quarks is not a straightforward task since 
quarks also (and mainly) undergo the strong interaction, in which calculations are not always 
accurate. Therefore, the SM prediction of many CKM observables suffers from uncertainties aris- 
ing from the use of approximations in QCD computations. These hard-to-quantify uncertainties 
do not belong to the well-defined framework of Gaussian statistics. 



"With more than three generations, there would be more CP violating phases. 

'The replacements p — > p = p(l — A^/2) and rj —y fj = rj{l — improves the accuracy of the UT apex. 



In order to avoid drawing wrong conclusions from the global CKM fit of such observables, 
one needs a robust and sufficiently conservative approach: we advocate the frequentist approach 
iZfitB, within the framework of the CKMfitter package, and describe its main features and results 
in the following. 



2 Statistical Framework and Fit Ingredients 

In the R&t frequentist statistical framework, theoretical uncertainties are treated without a 
priori knowledge except for the definition of rangesBty. In particular, no probability distribution 
functions are considered EJ. Practically, the fit likelihood is the product of an experimental part, 
^expixexp — 2;theo(ymod)) Comparing the measurements Xexp to the theoretical predictions xtheo 
(the latter depending on some model parameters 2/mod) 5 with a theoretical part, -Ctheo(?/QCD) f| 
which describes our knowledge on the QCD parameters ?/qcd £ ymod • A pseudo-x^ is built as 
= —2lnC and is minimized in the fit. 

The analysis is divided into three steps: first, the overall consistency between data and 
the theoretical framework (here, the SM) is tested by means of a Monte Carlo simulation. If 
the agreement between data and the SM is found to be acceptable, confidence levels (CL) in 
parameter subspaces are determined {e.g., in the p-f] plane). Finally, extensions of the SM can 
be tested and limits on new physics parameters are determined. 

Table |^ gives the input ingredientsof the fit. For \Vub\, the inclusive measurement of LEP0 
and exclusive measurement of CLEOEJ have been used, as well as the inclusive CLEO analysis 
using the B — > Xsj spectrumo. The input to the fit is the product of these three measurement 
likelihoods. The inclusive and exclusive LEP and CLEO \ Vcb\ measurements are averaged before 
entering the fit 111. For Am^ we translate the preferred value of the amplitude spectrum 
into a confidence level curve tZI. Finally, we use the world average of sin2P, dominated by the 
measurements of BABAR and Belle H. 



3 Standard Fit 

The consistency between data and the SM is gauged by the value of the overall minimum x^, from 
which a CL is inferred by means of a Monte Carlo simulation. Many sets of measurements are 
generated following the experimental likelihood distribution in which the theoretical predictions 
are computed for the best fitted values of the theoretical parameters. A fit is then performed on 
each of these sets, varying all the model parameters freely, and the CL of the SM fit is assigned 
by integrating the distribution thus obtained, up to the overall minimum x^- 

Figure ^ shows the distribution of the Monte Carlo generated (histogram), and the 
associated CL (smooth curve). The arrow shows the x^ of the standard fit, corresponding to a 
CL of 57%. The agreement between data and Standard Model is thus excellent, and we assign 
CL to various parameters of interest: Fig. |2| displays CL in the {p-fj) plane. For the individual 
constraints are shown the regions inside which CL > 5%, including the 4- fold ambiguity from the 
sin 2/? measurement. The latter is in perfect agreement with the SM prediction, demonstrating 
that the CKM mechanism is the dominant source of CP violation in fiavour changing processes 
at the electroweak scale. 



'^The four CKM parameters are part of these model parameters. 

"^A range of yqcD is obtained from QCD theoretical calculations; /^thco(j/QCD)=l for j/qcd inside the range, 
and when outside. 



Table 1: Input observables and parameters for the global CKM fit. When two errors are quoted, the first one is 
statistical and the second one systematic. When one error is quoted, it is statistical unless stated otherwise. 



Parameter 


Value ± Error(s) 


\Vud\ 






0.97394 ± 0.00089 


\Vus\ 






0.2200 ± 0.0025 


\Vub\ 


(LEP incl.)E 


_, (4.08 ± 0.61 ± 0.47) X 10"^ 


\Vub\ 


(CLEO incl.) 


11 (4.08 ± 0.56 ± 0.40) x lO^^ 


\Vub\ 


(CLEO excl.) 


B (3.25 ± 0.29 ± 0.55) X 10-3 


\Vcd\ 






0.224 ± 0.014 








A n/^n 1 A Aco 

0.959 ± 0.058 






/ A r\ /I 1 1 O 1 A A\ - , 1 A — 3 

(41). 4 ± 1.6 ± U.yj X lU 








(2.271 ± 0.017) X 10 ^ 


Anid 


le 




(0.496 ± 0.007) ps ^ 


Arus 


16 




See Text 


sin2/3wA0 


0.780 ± 0.077 


rric 






il.3±0.1syst) GeV 


mt(MS) 


(166.0 ± 5.0) GeV 


niK 






(493.677 ± 0.016) JVleV 


ArriK 




(3.4885 ± 0.0008) x 10 GeV 


mg^ 






(5.2794 ± 0.0005) GeV 


rriB, 






(5.3696 ± 0.0024) GeV 


mw 






(80.419 ± 0.056) GeV 


G p 






(1.16639 ±0.00001) X 10"^ GeV^^ 


Ik 






(159.8 ± 1.5) MeV 


Bk 






0.87 ±0.06 ±0.13 


Vcc 






1.38 ± 0.53sj^st 


Vet 






0.47 ± 0.04^y^t 


m 






0.574 ± 0.004,j,,t 


Vb{MS) 


0.55±0.01,j,,t 






(230 ± 28 ± 28) MeV 








1.16 ±0.03 ±0.05 



4 Beyond the Standard Fit 

4-1 Rare K Decays 

The detection of a second tt^uv event was announced early 2002 by the E787 Collabo- 

ration at Brookhaverta, leading to a branching ratio of (1.57lo'82) ^ 10"^'^ which is compatible 
with the SM prediction 113 given in Table |2|. A totaLof 5 to 10 SM events is expected to be 
detected at the successor experiment E949 at BNLE3, and thus this channel will become an 
important ingredient of the CKM fits. 

4-2 Charmless B Decays 

Both the a and 7 angles of the UT can in principle be extracted from the charmless two-body 
decays B — > vrvr and B Ktt , for which the data start to accumulate. Therefore, a lot of 
effort §§ has been put recently in the calculation of the charmless two-body hadronic matrix 
elements using the factorization concept. 
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Am^ (ps"') /(MC) 

Figure 1: Left figyfe: curves for various interpretations oLthe amplitude spectrum of Am^: the modified 
amplitude methodu (dotted line^the likelihood ratio methodtil (dash-dotted line), the expected x'^ curve for the 
current preferred value of Arris EZI (solid line). Right figure: Test of the goodness of the SM fit. The histogram 
(corresponding to the left scale) is the Monte Carlo-generated distribution (see text), which is integrated to 
compute a CL, as given by the smooth curve (right scale). The value of the current SM fit is indicated by the 

arrow, corresponding to a CL of 57%. 



Parameter 95% CL region Parameter 95% CL region 



A 


0.2221 ±0.0041 


\Vud\ 




0.97504 ± 0.00094 


A 


0.76 - 0.90 


\Vub\ [10- 




3.15 - 4.37 


P 


0.08 - 0.35 


[Kfcl [10- 




36.9 - 43.6 


V 


0.28 - 0.45 


\Vtd\ [10- 




6.3 - 9.1 


J [10-^ 


2.2 - 3.5 


\Vt,\ [10- 




36.4 - 43.0 


sin2Q; 


-0.81 - 0.43 


\Vtt\ 




0.99905 - 0.99932 


sin2/3 


0.64 - 0.84 


BR{Kl - 
BR(i4:+ 


TT°VV) [10-"] 


1.6 - 4.2 


a 


77° - 117° 


n+uD) [10-"] 


5.1 - 8.4 


P 


19.9° - 28.6° 


BR(B+ - 


T+Vr) [10-5] 


7.2 - 22.1 


7 = ^cp 


40° - 78° 


BR(J5+ - 


fi+iy,) [10-^ 


2.9 - 8.7 



Table 2: Fit results for the CKM parameters, the CKM matrix elements and branching ratios of some rare K 
and B meson decays, including the world average sin2/3 in the fit. Ranges are given for the quantities that are 

limited by systematic theoretical errors. 



One of these approaches, the QCD factorization approach (QCD FA) E3, has been imple- 
mented in CKMfitter. The results obtained from these calculations are labelled by "R&D", 
meaning that they are still subject to theoretical questions and that we do not intend to infer 
constraints on the CKM parameters from them. Rather, we intend to test the consistency of 
the QCD FA predictions in a global fit. 

In the QCD FA, the tree, strong penguin and electroweak penguin amplitudes are computed, 
the annihilation diagrams estimated, and the soft physics contribution in the spectator interac- 
tion are parameterized. The strong phases are found to be small, and thus QCD FA predicts 
small direct CP asymmetries. 

Global Fit of Branching Ratios and CP Asymmetries 

In a first stage, a global fit on QCD FA is performed using the recent charmless two-body 
branching ratios (BR), given in table ^, and time- integrated CP asymmetries (Aqp), given in 
Table ^ as measured by the BABAR , Belle and CLEO collaborations El. 

The numerous theoretical parameters are free to vary in the global fit according to the 
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Figure 2: Constraints in the p — f) plane: regions with more than 5% CL are shown for the main observables: 
Arrifj, Arria, ex, IKib/V'cbl- The global SM fit constraint lies at the intersection of all individual measurements. 
For both the global fit and the new world average sm2/3 constraint, two contour lines including points with more 
than 5% and 32% CL are drawn. Left plot: the sm2/3 constraint is not included in the global SM fit, but is simply 
overlaid: the measurement is in perfect agreement with the SM expectation. Right plot: the sin2/3 constraint is 
included in the global SM fit, and allows to pin down our knowledge on p ~ -q. 



Table 3: Branching ratios of the charmless two-body modes used in the global fit. The measurements from CLEO, 
BABAR and Belle collaborations E3 have been averaged into so-called World Average (WA) values. 



BR (xlO-*^) 


CLEO (9.1 fb-^ ) 


BABAR (55.6 fb'^ ) 


Belle (31.7 /ft"^ ) 


WA 


B° TT+TT- 
BO K+TT- 

S° ^ K+K- 
B+ ^ ^+^0 
B+ ^ K+ttO 


4.3li;^±0.5 

n.2t]:l ± 1.2 

< 1.9(90%) 
5.6+2 3 ± 1-7 

11 R+3.0-1-1.4 
ii.D_2.7„i.3 

18.24-1 ± 1.6 

i4.6i:f4.3 


5.4 ±0.7 ±0.4 
17.8 ± 1.1 ±0.5 
< 2.5(90%) 

5.1 ±2.0 ±0.8 
10.8 ±2.1 ±1.0 
18.2 ±3.3 ±2.2 

8.2 ±3.1 ±1.2 


5.1 ± 1.1 ±0.4 
21.8 ± 1.8 ± 1.5 

< 2.7(90%) 
7.0 ±2.2 ±0.8 
12.5 ±2.4 ±1.2 
18.8 ±3.0 ±1.5 
7.7 ± 3.2 ± 1.6 


5.17 ±0.62 
18.5 ± 1.0 

5.9 ± 1.4 

18.5^2 2 

9.0 ±2.2 



predications of Ref. E3. In the left plot of Fig. ^, the QCD FA constraints from the charmless 
two-body BR and Aqp are compared to the SM constraints coming from all the measurements 
of Table |l|. 

One concludes that the QCD FA can accommodate the branching ratios and CP asymmetries 
of the i? — >■ TTTT and B — >■ K-k decays, and that the preferred value of the angle 7 is around 80°. 
The constraints in the p-f] plane are in agreement with those obtained in the standard fit. 

Ctttt and 5^7r Fits using QCD FA 

The time-dependent CP asymmetry of B^ — > tt^tt^ rate reads 



acp{t) = S^TrSin{AmdAt) - C^^cos(AmrfAt) , (1) 




-1 -0.5 0.5 1 



P p 

Figure 3: Left Plot : CL in the (p — 77) plane obtained from a fit using branching ratios and CP asymmetries 
from Table ^, in the framework of the QCD FA approach. Right Plot: CL in the (p — fj) plane obtained from a 
fit using C-rnr and St„v from BABAR where S and the ratio |P/T| are predicted by the QCD FA approach. The 
light, medium and dark shade zones include points with more than 5%, 32%, 90% CL. 

Also shown is the 95% CL zone of the SM fit. 



with At, the time difference between the CP and the tag B decays, and with the coefficients 
given by 

_ 2ImA7r7r An— ^ ~ l-^^^l^ /'o^ 

i ~r 1^71-77 1 i -r IAttt^I 

The CP parameter Atttt is given by 

_ 2^0 A(B' ^ 77-^77-) _ 1 - {Rt/ Ru)\P/T\e'^^--'^ _ 

^(^0^^+^-) l-(i?t/i?„)|P/r|e+^(°+^) " ' ' ^ ' 

where the phase e~^*^ arises from B^'B^ mixing, Rt(u) ^-^e the sides of the UtI, and T, P are 
the tree and penguin amphtudes. In the presence of penguin contributions, the UT angle a 
enters the expression of Atttt with a correction term which depends, in particular, on the relative 
strong phase 5 = arg(PT*) between the penguin and the tree amplitudes. As a consequence, a 
measurement of the parameters S-j^t^ and Ctttt cannot be translated into sin2a without information 
on the relative strong phase 5 and the ratio \P/T\, e.g., as computed by the QCD FA approach. 

The right plot in Fig. ^ and both plots in Fig. ^ show CL in the p — fj plane using the 
Ctttt and /Stttt measurements of BABAR and Belle respectively. While the BABAR measurement, 
interpreted within the QCD FA approach, is in good agreement with the SM expectations, the 
Belle measurement is in strong disagreement. 

5 Beyond the Standard Model: the example of j4sl 

Two classes of models beyond the SM have been investigated @ using the semi-leptonic asym- 
metry which measures CP violation in iJ^-mixing: 

, ^ r(i30(t) ^ l+X) - T{B\t) ^ l-X) _ Fia 
T{B^{t) ^l+X)+T{B^{t) ^l-X) ™Mi2' 



Table 4: Time-integrated direct CP asymmetries A in B ^ Kn and CP parameters Cttt and Stt-k of the 

time-dependent analysis in _B" tt^tt". 



Observable 


Value 


ACP(B^ K+TT-) 

Acp{B+ ^ K+tt'^) 
Acp{B+ K°7r+) 


-0.05 ± 0.05 (WA) 
-0.09 ±0.12 (WA) 
0.18 ±0.10 (WA) 


<? 
c 


-0.01 ± 0.38 (BABAR) 
-1.2ll°;|J (Belle) 

-0.02 ± 0.30 (BABAR) 
-0-94 ta27 (Belle) 




Figure 4: Left Plot : CL in the p — f) plane obtained from a fit using Crvr and 5,rir from Belle, where the relative 
strong phase 5 and the ratio |-P/T| are predicted by the QCD FA approach. The light, medium and dark shaded 
zones include points with more than 5%, 32%, 90% CL. Also shown is the 95% CL zone of the SM fit. Right Plot: 

zoom of the CL between and 10%. 



The present world-average cZl reads ^SL = (0.2 ± 1.4) x 10~^. This parameter is predicted to 
be small (of the order of 10"'^) within the SM, but can be easily enhanced by new physics 
contributions. 

In a class of models beyond the SM in which the CKM matrix is still a 3 x 3 unitary matrix 
and in which tree processes are dominated by the SM contributions (thus = r^2^, and only 
Mi2 is affected by new physics), M12 can be parameterized by 

M12 = r^^^Mf^, (4) 

where and 6d are a modulus and a phase respectively. Constraints from sin 2/3, lS.Md and Asl 
are thus also modified. 

The two left plots of Fig. ^ show CL in the (r^ — 20^^) plane using all present constraints of 
the Bd meson system, including (middle plot) or not (left plot) the constraint from Asl: one 
sees that the ^SL constraint disfavours the small and large sm.29d values. 

A particular subset of the previous class of models have no new phase in flavour changing 
processes: these are the Minimal Favour Violating (MFV) models, as the minimal MSSM, the 



two-Higgs doublets model, etc. In these models, one has: 



where So{xt) is the Inami-Lim functionEl depending on the top quark mass. Ftt enters AM^, 
AMs, ex, sin 2/3 and Asl- The right plot of Fig. ^ shows CL in the {p — fj) plane within the 
MFV framework, where negative values of fj are no longer excluded ^. As can be seen through 




Figure 5: Left Plot : CL in the — 29d plane using the present standard constraints of the Bd system. The 
light, medium and dark shaded zones corresponds to CL > 10, 32 and 90% respectively. Middle Plot : CL in the 
— 26d plane using the standard constraints of the Bd system as well as the Asl measurement. Right Plot: CL 
in the p — plane in the framework of the MFV models. Negative values of f] are allowed. 



the example of ^SL, a global CKM analysis is not only the place to test the Standard Model 
CKM paradigm, but also to investigate New Physics models, on which one can start to put 
meaningful constraints. When measurements reach a better precision, one will be able to rule 
out models, while favouring others. 



6 Conclusion 

The successful operation of the B factories at SLAC and KEK has provided a flood of new 
results that can be used to obtain information on the CP- violating phase of the CKM matrix. 
The extraordinary agreement of the world average of sin2/3 with the indirect determination 
represents a new triumph for the Standard Model. It establishes the KM mechanism as the 
dominant source of CP violation in flavour changing processes at the electroweak scale. While 
the inclusion of sin2/3 into the global CKM fit is a theoretically straightforward task, exploiting 
rare B decays requires some further theoretical understanding in order to lead to meaningful 
constraints at present. 
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